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ABSTRACT 

We explore the impact of Fermi-like acceleration of Lyman-alpha (Lya) photons across 
shock fronts on the observed Lya spectral line shape. We first confirm the result of 
Neufeld & McKee (1988) that this mechanism gives rise to extended blue wings which 
may have been observed in some radio galaxies. Our Monte-Carlo radiative transfer 
calculations further show that in a minor modification of the shell-model, in which 
we add an additional static shell of hydrogen, this process can naturally explain the 
small blue bumps observed in a subset of Lya emitting galaxies, which have been 
difficult to explain with conventional shell-models. Blue bumps can be produced with 
an additional column density of static hydrogen as small as ^ 

typically occur at roughly the outflow velocity of the shell. In our model the spectra of 
so-called ‘blue-bump objects’ might reflect an evolutionary stage in which the outflows 
regulating the escape of Lya photons are still engulfed within a static interstellar 
medium. 

Key words: line:proflles — radiative transfer — shock waves — methods:numerical 
— galaxies: star formation — ISM: jets and outflows 


1 INTRODUCTION 


Observations indicate that the escape of Lya photons is fa¬ 
cilitated enormously by the presence of outflowing interstel¬ 
lar gas (e.g. Kunth et al.|1998 ). Scattering off these outflows 


effectively Doppler shifts Lya photons out of resonance into 
the red wing of the absorption line profile. In models of this 
process the outflow is often represented with a geometrically 
thin shell of gas (e.g. Verhamme et al.||2006 Schaerer et al. 
[2011 ). In spite of its simplicity, this so-called ‘shell-model’ 


has been very successful at reproducing observed spectra 
line shapes ( Verhamme et aL]|20Q8 Hashimoto et al.|2015 ). 

However, the shell-model has recently been shown to 
have difficulties, especially in reproducing the strength of 


‘blue bumps’ in a subset of observed spectra (Kulas et al. 


2012 Chonis et al. 2013 Hashimoto et al. 2015). Adams 


et al. (2009) studied what may be viewed as an extreme ex¬ 


ample of blueshifted Lya emission in which the Lya line of 
spatially extended Lya emission around a radio galaxy as a 
whole is blueshifted relative to the observed 21-cm absorp- 
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tion, with emission extending beyond 1500 km s ^ into the 
blue wing of the line profile. 


Blueshifted Lya radiation arises naturally when the 


photons scatter through optically thick inflowing gas (Zheng 

& Miralda-Escude 

2002 

Dijkstra et al. 

2006a|b). Indeed, 

Adams et al. (2009 

) show that their data can be reproduced 


large mass of cold neutral gas inside a massive dark matter 
halo (M ^lO^^M©) though, is at odds with expectations 
from theory, which predicts that the gas should be predom¬ 


inantly accreted in the hot-mode (e.g. Keres et al. 
Fig. 7 of Dekel Birnboim||2006 ). 


2005 


An alternative process which gives rise to blueshifted 
Lya emission is described by Neufeld & McKee (1988; hence¬ 
forth NM88), who show that scattering across a shock front 
can give rise to such blueshift. NM88 presented analytical 
calculations, and therefore were forced to adopt simplifying 
assumptions, namely that the shock-crossing acceleration 
process is modelled by Lya photons bouncing between two 
partially transparent, frequency-preserving mirrors. In real¬ 
ity, however, the frequency of a Lya photon is not preserved 
at each ‘reflection’. Instead, the frequency of the photon dif¬ 
fuses through frequency-space as the photon is scattered by 
the material on either side of the shock front, in turn chang¬ 
ing the optical depth of the gas to the photon. In addition. 
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2 Chung et al. 


they assume that photons cross the shock front isotropically. 
Finally, the analysis presented by NM88 was restricted to a 
simplified geometry of two adjacent semi-infinite slabs. 

In spite of an increasing prevalence of Lya radiative 


Cantalupo et al.||20051 Dijkstra et al.||2006al 

Tasitsiomi, A. 

2006| IVerhamme et al. 2006| [Laursen et al. 

|2009| Yajima 


et al. |2012), they have so far not been used to investigate 


this Fermi-like acceleration mechanism (while acknowledg¬ 
ing the inaccuracy of doing so, we henceforth refer to this 
as ‘Fermi acceleration’ for brevity). The goals of this pa¬ 
per are two-fold: (z) numerically study Fermi acceleration of 
Lya photons across a shock front without the simplifying 
assumptions that were required in the analytical treatment, 
and (a) investigate whether this mechanism can help to re¬ 
produce observed blue bumps in spectra that have been dif¬ 
ficult to explain with conventional shell-models. The outline 
of this paper is as follows: § [^describes the basic acceleration 
mechanism, § |3.1| details the simulations which we use to in¬ 
vestigate the physical mechanism, § |3.2| presents the output 
of our numerical experiments, discusses the applications 
of our findings, and in §[^we summarise the conclusions. 


boost will impart a blueshift in the local gas frame. Par¬ 
tially coherent scattering off atoms with thermal velocity 
Vth mostly preserves the blueshift of this photon: for a pho¬ 
ton with frequency Xs ~ Vs/vth each scattering event pushes 
back the photon to the line centre by an average amount 
— l/|xs| ( Osterbrock||1962| . When the photon scatters back 
across the shock front before this ‘restoring force’ has re¬ 
turned the photon to line center, the blueshift of the photon 
increases with each shell crossing. This is what NM88 re¬ 
ferred to as ‘Fermi acceleration’. 

In order for Fermi acceleration of a Lya photon to 
actually occur we require it to scatter after crossing the 
shock front. The optical depth to a Lya photon that is 
blueshifted to Vb through a slab of gas with column density 
Nm equals r(vh) ~ 0.6(A^hi/ 10^^ cm“^)(u6/200 km 
The condition r{vb) ^ 1 thus translates to N ^ 

1.6 X 10^^(u6/200kms“^)^ cm“^. Inversely, a photon can 
be Fermi accelerated to a maximum blueshift of Vb ~ 
155(A^hi/ 10^^ cm“^) km s“^, in the frame of the seattering 
medium. After scattering, the Lya photon can get an addi¬ 
tional Doppler shift ~ Ug depending on whether the scatter¬ 
ing medium is moving or not. 


2 BLUESHIFTING MECHANISM 

Lya photons are resonantly scattered by HI, making the 
Lya radiative transfer complex. The propagation of Lya 
photons is affected by both bulk gas motion and, in detail, 
microscopic motion of individual hydrogen atoms. The re¬ 
sult is that Lya photons undergo a random-walk like motion 
in both physical and frequency-space (see Dijkstra|2014| for 
a detailed review). 

As Lya photons diffuse away from the resonance fre¬ 
quency, their mean free path - and hence their escape prob¬ 
ability - increases. As a result, the spectrum of Lya photons 
emerging from static optically thick media consists of two 
peaks, which are distributed symmetrically around the res¬ 
onance frequency (as frequency diffusion can occur to lower 
and higher energies with equal probability^ . When the scat¬ 
tering medium is contracting [expanding] however, the fre¬ 
quency diffusion preferentially occurs towards higher [lower] 
frequencies. As a result, the spectrum of Lya photons emerg¬ 
ing from contracting optically thick media is blueshifted. An 
alternative way to see this is that the converging flow of 
contracting gas is doing work on the Lya photons as they 
scatter outwards, which increases their mean energy (Zheng 
fc Miralda-Escud^|2002 Dijkstra et al.||2Q06a|b ). 


The Fermi acceleration mechanism described in NM88 
invokes similar converging flows of gas. The mechanism is 
illustrated with a simplified geometry in Figure which 
shows two semi-inhnite adjacent slabs of neutral gas. The 
slab on the left is moving to the right with velocity Vs into a 
stationary slab. If Vs exceeds the sound speed in the left slab, 
then the two slabs are separated by a shock front. When a 
Lya source is in the vicinity of this shock front we expect 
some Lya photons to diffuse through space and cross the 
front. When a photon traverses the shock front, the Doppler 


^ For low gas temperatures (T = 10 K) the energy deposited 
in the recoiling scattering atom becomes important, and the red 
peak is enhanced. 


3 FERMI ACCELERATION IN 
MONTE-CARLO SIMULATIONS 

In this Section we describe the setup and results of the nu¬ 
merical simulations which we use to study the Fermi accel¬ 
eration mechanism briefly summarised in § [^ 


3.1 Simulation Setup 

We use modified versions of the SLAF 3D Monte-Carlo Lya 
radiative transfer code (Chung et al. 2015, in prep), as well 
as the code described in Dijkstra et al. (2006a). The vanilla 


version of SLAF is a grid-based code, able to handle arbi¬ 
trary gas and source distributions within a finite volume 
of 3D space. For this work we have modified it, removing 
the dependence on an underlying grid-structure, so that the 
code can handle the specific setup described in NM88. 

We choose to focus on the geometry discussed in NM88 
to facilitate a straightforward comparison. Explicitly, there 
are two distinct regions of gas, represented by slabs of infi¬ 
nite width (y) and height (z). The slab on the left side repre¬ 
sents outflowing gas, and has a bulk velocity of Vs = 400 km 
s“^ from left to right. The HI column density across the out¬ 
flowing slab is Al = 1 X 10^°cm“^. The slab on the right side 
represents static gas with an HI column density Nc = 10^^ 
cm“^ ^ Ni. A shock front exists at the interface of the 
two slabs. Lya photons are emitted at line centre in the rest 
frame of the outflowing slab on the left, which represents 
Lya emission powered by star formation triggered by the 
passing shockwave (see NM88). We ignore dust throughout 
the analysis in § |3.2| which simplifies the interpretation of 
our results and does not affect our main results. We discuss 
the impact of this assumption separately when we discuss 
applications of our results in § |4.1| and § |4.2[ Figureshows 
the described experimental setup. 

Figure 2 was produced from a run tracking 5 x 10^ cross¬ 
ing events, while Figures 3 and 4 were produced from runs 
of 5 X 10^ photons. In all cases the results were checked for 
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Figure 1. Schematic showing the setup used for studying the 
basic acceleration mechanism (see § and 1^. The Lya source 
is denoted by the cross in the left slab. The left slab has HI column 
density A/^i, and the right slab has HI column density Nc- The gas 
in the left slab moves with velocity Vg, photons cross the shock 
front with velocity Vp^otoni ^-nd the angle between Vg and Vp^oton 
is labelled 6. is then the projection of Vp^oton onto the unit 
direction vector of the shock front, Vg. 


convergence, and found to already be converged with much 
lower photon counts. 


3.2 Results 

3.2.1 Shock Crossing Statistics 

Table shows the fraction of photons escaping the simula¬ 
tion volume, /n, as a function of the total number of shock 
crossings, n. The second and third columns indicate the fn 
found in our simulations and the analytic estimate, respec¬ 
tively. The following paragraphs explain how we obtained 
these analytic estimates. 

Photons emitted in the centre of the outflowing slab 
are equally likely to leave this slab on the left or right hand 
side. We therefore expect /o = 1/2 of all photons not to cross 
the shock front, which is in agreement with the simulated 
fraction. 

Photons that enter the static slab for the first time are 
Doppler boosted by an average amount Vb = 2vs/3 (see 
§ 3.2.2). The probability of these photons being transmit¬ 
ted all the way through the slab is given by T = 4/(Sri), 
where Ti denotes the optical depth of the entire static slab 
to photons that enter at frequency Xi (in the frame of the 
slab, see |Neufeld||1990| note that this transmission proba¬ 
bility only applies when the optical depth of the slab to 
incoming photons is ^ 1). Substituting numbers yields 
T O.Obli’^y/A^Hi? where A^hi denotes the HI column den¬ 
sity of the slab in units of 10^° cm“^, and Vb? denotes Vb 
in units of 100 km s“^ (notation adopted from NM88). We 
then get /i = (1 - /o) X T 0.018, which compares reason¬ 
ably to the fraction fi = 0.023 we found in our simulation 


Table 1. Fraction of photons undergoing n successive shock cross¬ 
ings before exiting the simulation (middle column). The right 
column refers to the same fraction calculated with an analytic 
method. 


total # shock crossings (n) 

fraction (fn) 

analytic estimate 

0 

0.50820 

0.5 

I 

0.02299 

0.018 

2 

0.26215 

0.28 

3 

0.04578 

0.068 

4 

0.12525 


5 

0.01283 


6 

0.01944 


7 

0.00149 


8 

0.00168 


9 

0.00008 


10 

0.00011 



(exact agreement is not expected as photons cross the slab 
over a range of frequencies). 

Photons that cross the front twice are back in the low 
column density outflowing slab, and appear Doppler shifted 
by an average amount Vb ~ 4vs/S in the slab frame (see 
§ 3.2.2). The slab will appear optically thin to these pho¬ 
tons, and we estimate the transmission probability from 
T ~ exp(-2Ti) + 0.5[1 - exp(-2Ti)]. Here, n = Nia^Vb) 
denotes the optical depth through the outflowing shell for 
the incoming photons, where acx(vb) is the Lya absorption 
cross-section at Vb. The factor of 2 in the exponent accounts 
for the fact that the photons enter the slab under an an¬ 
gle 0 with probability P(cosO) oc cos^ (see § 3.2.2). The 
transmission probability is thus the sum of the probability 
(z) that photons are transmitted directly through the slab 
(the exp(—2Ti) term), and (ii) that photons that do scat¬ 
ter eventually escape the slab without crossing the shock 
again. For the latter, we assume that photons are equally 
likely to escape on the left and right hand side after scat¬ 
tering, which seems reasonable given that the slab appears 
optically thin on average. Substituting numbers we obtain 
/2 = (1 — /o — /i)T ~ 0.28, which agrees with the simulation 
result. 

With T = 4/(Sri) and Vb ~ 2vs, we estimate fs = 
(1 — fo — fi — / 2 ) X T ~ 0.068, 48.5% larger than what we 
obtain from the simulation. We attribute this discrepancy 
to the approximation T = 4/(Sri) breaking down at larger 
shifts from line centre. 

The previous analysis allows us to understand quanti¬ 
tatively the simulation results. More photons escape after 
undergoing an even number of shock crossings than odd. 
This is simply a reflection of the fact that in our experimen¬ 
tal setup we have specified that Nc^ Ni. 

As we pointed out in §[^ it is theoretically possible that 
a Lya photon that has been Fermi accelerated into the line 
wing can subsequently diffuse back into the core of the line 
prior to escaping or crossing the shock front again. When 
this happens, any memory of previous Fermi acceleration is 
erased. However, our numerical simulations indicate that in 
practice this almost never occurs, and therefore we can iden¬ 
tify odd [even] numbered contiguous shock crossings with 
photons which exit on the right [left] side of the simulation. 
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Figure 2. Distribution of the projection of the shock-crossing 
angle for various values of 1. For each value of I the plotted distri¬ 
bution includes all photons that make an crossing regardless 
of what their final crossing count, n, is. 


3.2.2 Velocity Shift vs Shock Crossing Number 


We denote the angle at which a photon crosses the shock 
front with 0 (see Fig [^. For an isotropic distribution of 
shock crossing directions we have P(/i) = 1, where = cos 0. 
The average Doppler shift experienced by a photon as it 
crosses the shock front is Vb = Vs fg dfi = Vs/2. After 

I crossings the photon experiences an average blueshift of 
lvs/2 (NM88). 


With our code we are able to track individual photons 
as they cross the shock front, which allows us to directly 
measure Figure]^ shows the number count, for 

photons that cross the time for various 1. P(/i) only dif¬ 
fers from N(fi) by a normalisation factor. Here, I refers to 
the current (rather than final) contiguous crossing count of 
a particular photon as it crosses the shock front. A photon 
which finally exits the simulation after n crossings is rep¬ 
resented n times (once for each crossing). Figure shows 
clearly that the distribution of photon crossing projections 
is not isotropicj^ Instead, we find that oc //. This dis¬ 


tribution has been found in previous analyses (e.g. Ahn et 
[ar]|200l[ [Garavito-Camargo et al. 2014| , and its origin is 
discussed in Appendix |A| For this di strib ution we expect 
Vs fg djj. — 2vs/3 (see also Fig. Al), and hence that 

a photon experiences a Doppler boost 2lVs/3 after I shock 
crossing^ 

Figure shows the velocity shift Vb of exiting photons 
as a function of the total number of shell crossings n, in 
the frame of the gas into which the photon is crossing 
(i.e. for odd [even] numbered shock crossings the velocity 
shift is given in the static [outflowing] frame). Results 
from our Monte-Carlo simulations are represented by the 
data points. The red dashed line shows the analytic result 


^ In Figure p, > 0 [/i < 0] for photons that cross the shock 
front from left-to-right [right-to-left]. In Figure]^ only photons 
on an odd [even] shock crossing contribute to N(/j,) for positive 
[negative] values of g. 

^ Note that if a photon crosses the shock front at an angle 0, 
then the total hydrogen column density to the edge of the slab 
is Nm/g^ where Ahi = Ni [A^hi = Nc] for the slab on the left 
[right] in Figure^ The angle-averaged column to the edge of the 
slab is given by (A^hi) = Jq dg P(g)Niii/g = 2Ahi, as used in 
§ 3.2.1. 


Figure 3. Mean velocity offset for photons exiting with n shock 
crossings as found in our simulations (points), as measured in 
the frame of the gas. The error bars show the standard devia¬ 
tion within the distribution. The dashed red line shows the ^Ivg- 
relation which applies to isotropic shock crossings, whereas the 
green line refers to the |/4;s-relation which applies to anisotropic 
shock-crossings of the form P{g) oc g (see text). 


under the assumption that photons cross the shock front 
isotropically (as in Neufeld McKee|1988 ), while the green 
solid line shows the analytic result for photons crossing 
the shock front according to P(/i) oc /x. For n ^ 3 the 
simulation results follow the -relation (where here 
n is a good proxy for /), after which it approaches the 
isotropic shock-crossing case. The reason for this transition 
is that as the photons get blueshifted further into the wing 
of the line profile their mean free path increases. As we 
discuss in detail in Appendix [A] we expect a transition to 
isotropic shock crossing when the mean free path becomes 
comparable to the thickness of the slab. 


Finally, Figure ^ shows the spectrum outcoming from 
the slab, as measured in the lab frame. The result can be 
easily understood in terms of our previous analysis. More 
specifically, the total spectrum (black histogram), can be 
explained as follows. 


• The peak at Vb ~ 300 km s~^ is composed primar¬ 
ily of the /o = 50% of all photons that did not cross 
the shock front. These photons diffused outward of the left 
slab, and would have emerged with a characteristic dou¬ 
ble peaked emission line profile centered around the reso¬ 


nance frequency (e.g. Adams 1972 Harrington 1973 Neufeld 


1990). However, when we Doppler boost these photons back 


into the lab frame this double peak is diluted by the fact 
that the photons escaping the slab do so at various angles. 
The Doppler boost is dependent on the exit angle, with an 
average Doppler boost of {/j)vs ~ 270 km s“^. The grey 
crosshatched histogram shows the spectrum of these pho¬ 
tons. 

• The small fraction of photons that escape after a sin¬ 
gle shock crossing, fi ~ 2%, is indicated by the purple 
hatch-filled histogram and escapes with a mean blueshift 
of ~ 2vs/3 ~ 300 km s~^ (also see Fig. i these photons es¬ 
cape from the static shell, and no additional Doppler boost 
into the lab frame is required). 

• The photons that escape after two shock crossings 
(green dotted line histogram) account for /2 ~ 26% of the 
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Figure 4. Spectrum of photons exiting the slabs with a single 
source, lab frame, n is the number of shock crossings that the 
photon undergoes prior to exiting the simulation (see § |3.2.1| ). 

total. Figure shows that these photons have accumulated 
a blueshift of ~ 4us/3 ~ 530 km s“^ in the frame of the out¬ 
flowing gas. A Doppler boost back into the lab-frame trans¬ 
forms these photons back to a blueshift of ~ 2us/3 ~ 300 
km s“^. 

A similar reasoning applies to photons that escape af¬ 
ter a larger number of crossings, and shows that it is easy to 
understand the shape of the total spectrum: the redshifted 
peak (grey crosshatched histogram) at Vh ~ 2us/3 is com¬ 
posed of photons that did not cross the shock front. The 
first peak blueward {vh < 0) of the Lya resonance is at 
Vh ~ —2vs/S and consists of photons that crossed the shock 
front 1 — 2 times. Similarly, the second blueward peak at 
Vh ~ —4us/3 consists of photons that crossed the shock 3-4 
times, etc. The prevalence of these peaks depends on a num¬ 
ber of factors, including Vs, iVi, Nc, and also the distribution 
of Lya sources relative to the HI gas. 

Finally we note that our results here are in excel¬ 
lent agreement with the analytic estimate of the maximum 
blueshift detailed in §[ 2 ] Setting N = 10^^ cm“^ in the ana¬ 
lytic formula gives Vb ~ 1600 km s“^, which lies in the blue 
tail of Figure 


4 APPLICATIONS 

We now explore, by means of our Monte-Carlo codes and 
assuming some simple models, whether the Fermi accelera¬ 
tion mechanism can provide a viable explanation for some 
observed systems. 


4.1 Extended Blue Wings in Radio Galaxies 

The analysis by NM88 was motivated by observations of 
3C326.1, a radio source with a Lya spectral line profile with 
emission extending far into the blue wing of the line pro¬ 
file ( Djorgovski||1988 Strauss et al.|p'987 ). NM88 modelled 
3C326.1 with a shock front that was propagating into a col¬ 
lection of dense clumps. The shock passage triggers star for¬ 
mation and Lya emission inside the post-shock gas. The sub¬ 
sequent radiative processes in each clump proceed just as in 
the slab models we discussed previously. The details of the 


emerging spectrum depend on the distribution and covering 
factor of the clumps, but for simplicity and to facilitate a 
direct comparison to NM88 we follow their prescription and 
assume that the emerging spectrum is a superposition of 
the spectra of individual clumps. NM88 further show that 
dust does not affect the emerging spectra, provided that 
the gas-to-dust ratio ^dust ;^0.016, where ^dust denotes the 
dust-to-gas ratio relative to the local interstellar value. They 
argue that this limit is acceptable because the dust grains 
would be destroyed by the passage of the shock front. We 
finally assume all clumps - and therefore the radiative trans¬ 
fer processes inside them - to be identical, then we can sim¬ 
ply adopt the spectrum from Figure In reality, from our 
findings in § 3.2.2, we expect that a distribution of clump 
properties would give rise to a superposition of spectra with 
peaks in different locations, and therefore that the emission 
in the blue wing of the line profile would have less prominent 
feature^] 

Putting these reservations aside, we observe that qual¬ 
itatively our simulated spectrum shows similarities to that 
predicted by NM88. We should note though that in Figure]^ 
we omit the direct sources at = 0, which NM88 put in by 
hand as a delta function at = 0. 

We note that the two largest peaks immediately red- 
ward and blueward of line centre cannot be attributed to 
photons exiting without crossing the shock front as in NM88 
(i.e. the standard double-peak profile produced by a static, 
optically thick medium). Instead, our analysis shows clearly 
that these peaks are primarily caused by photons which have 
undergone 0 and {1, 2} shock crossings respectively. The 
smaller, secondary blueward peak is primarily composed of 
photons which have undergone 3 or 4 shock crossings. The 
primary differences between the spectra predicted by us and 
by NM88 are in the asymmetry of the two main peaks, and 
the ratio of the peaks. As we have previously mentioned, the 
source distribution and physical parameters of the clumps 
affect the out coming spectrum, and it is likely that tuning 
of either or both of these will alleviate this discrepancy. The 
two blue peaks in our model appear to be shifted redward 
compared to the model predictions from NM88. However, 
a comparison of the observed spectrum in NM88 and our 
model shows that the positions of our blue peaks are also a 
good fit to the data. In fact, if the amplitude of the peaks is 
ignored (this is dependent on the exact source distribution), 
the position of the peaks in our model is arguably a better 
fit to the data than the NM88 model. 

We now consider the observed Lya spectra of z ~ 3.4 
radio galaxy B2 0902+34 (Adams et al. 2009). There is 
broad agreement in the shape of the observed and pre¬ 
dicted spectra, with the broad wing extending well beyond 
~ 1000 km s“^, and some spectra displaying prominent 
peaks blueshifted by ~ 1000 km s“^. The spectra from B2 
0902+34 vary widely depending on where exactly on the 
galaxy the fiber is placed. It is however true to say that for 


Even for fixed clump properties, the emerging Lya spectrum 
depends on the placement of Lya sources. In particular, chang¬ 
ing only the source position leaves the peak positions unchanged, 
but modulates their relative amplitudes. We have explicitly ver¬ 
ified, however, that our spectrum barely changed if we assumed 
a uniform distribution of Lya sources throughout the outfiowing 
slab. 
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the majority of fibers the most prominent peak has been 
observed to lie either at the systemic velocity, which was 
determined from the observed 21-cm absorption signature, 
or slightly blueward of it. This would appear to be in ten¬ 
sion with our model here which, following NM88, predicts 
the most prominent peak to appear at I’b = 0 (not shown in 
Figure]^ as a result of direct ‘blister’ sources on the edge of 
the clump. Recalling that the amplitude ratios of the peaks 
is determined by the source distribution, we speculate that 
this tension could be relieved if the most prominent peak 
in the B2 0902+34 observations is instead identified as the 
first blueshifted peak in our model, and the absence of a 
pronounced peak at systemic velocity attributed to a lower 
prevalence of ‘blister’ sources in this system. 



4.2 Blue Bump Spectra 

As we mentioned in § ^ observed Lya spectral line profiles 
can often be reproduced surprisingly well with shell-models. 
In these models, a central Lya source is surrounded by a geo¬ 
metrically thin shell of gas. The shell-models contain two pa¬ 
rameters describing the Lya source: (z) the assumed FWHM 
of the Lya line prior to scattering, and (ii) the strength of 
the Lya emission line, which is quantified by the equivalent 
width (EW). The shell itself is described by four additional 
parameters: (z) the HI column density of the shell (zz) 

its outflow velocity z;sheii, (Hi) its dust content ra, and (zz;) 
its velocity dispersion b. 

Some recent analyses have pointed out that shell-models 
have difficulties explaining blue bumps in observed spectra 


(e.g. Kulas et al.|2012 Chonis et al.|20l'3 Hashimoto et al. 


2015 ). Nice examples can be seen e.g. in Figure 7 of Chonis et 
al. ( |2Q13 ). More recently Hashimoto et al. (2015) pointed out 
that the shell-models, when applied to blue-bump objects, 
consistently require the intrinsic Lya line to have a FWHM 
which is too large (an excessively high FWHM was only 
required for blue-bump objects). 

Because Fermi acceleration naturally gives rise to 
blueshifted emission, we investigate a minor modification 
of the shell-model in which we embed the outflowing shell 
within a static gas cloud. This modification can be inter¬ 
preted as the situation in which the outflow is still prop¬ 
agating into the static interstellar medium (possibly prior 
to breaking out of it). This configuration now includes a 
shock front as in our previous analysis. We specifically study 
a shell-m odel with parameters bas ed loosely on those in¬ 
ferred by 


Hashimoto et al. 


(|2015|: 


= 10 ' 


^sheii = 200 km s“^, and FWHM=200 km s“ . We fur¬ 
ther assume ra = 0, EW= oo (i.e. pure Lya emission), and 
b = 12.9 km s“^ (corresponding to gas at lO^K), but note 
that these assumptions do not affect our results at al|^ Ei- 
nally, the key new model ingredient is the addition of a static 
shell of gas adjacent to the outflowing shell. This shell is 


^ For completeness we have presented our numerical results which 
include dust in Appendix These results show clearly that dust 
barely affects the blue bumps. This result is not surprising: our 
modification only adds a small amount of additional hydrogen, 
and this additional gas triggers Fermi-acceleration into the (blue) 
wings of the Lya line profile, where Lya photons escape more 
easily. 


Figure 5. Predicted Lya spectra emerging from a Lya source 
surrounded by a shell of HI gas outflowing at 200 km s“^ (the 
shell-model, see text) embedded within a static neutral shell with 
a column density . The solid black histogram shows the spec¬ 
trum for =0 (i.e. the standard shell-model), while the red 

dashed and blue dotted line shows the spectrum for = 10^® 

cm“^ and = 5 x 10^® cm respectively. The plot demon¬ 

strates that adding even a small amount of hydrogen (1/10^^ of 
that in the shell) triggers the on-set of Fermi acceleration, which 
gives rise to a blue bump. 


characterised primarily by its column density N^i^. We as¬ 
sume that it has no dust and that it has the same tempera¬ 
ture as the outflowing shell. 

The results of this analysis for = 0 (i.e. the orig¬ 
inal shell-model; solid black histogram), = 10^® cm“^ 

(red dashed line), and = 5 x 10^® cm“^ (blue dotted 

line) are shown in Figure]^ where it is clear that a small ad¬ 
ditional column of hydrogen dramatically 

affects the spectrum blueward of the systemic velocity. This 
large change of the spectrum can be easily understood, as 
the outflowing shell with — 10^^ cm“^ directly trans¬ 

mits a significant fraction of Lya photons. However, the sur¬ 
rounding static shell is optically thick to photons emitted 
near line centre, because these photons still appear close to 
the centre of the line in the frame of this gas. The static 
shell therefore effectively reflects back photons into the out¬ 
flowing shelj^ The reflected photons appear blueshifted by 
~ Vshell in the frame of the outflowing shell, where their 
newly acquired large blueshift makes them escape efficiently. 
For example, photons that are scattered by 90° after be¬ 
ing reflected back into the outflowing shell escape with a 
blueshift of ~ z;^, which is indeed where the new peaks in the 
spectrum emerge. This process is depicted schematically in 
Figure The analytic estimate for the maximum blueshift 
in § [^ provides a decent estimate of the maximum blueshift 
Vb ~ 160kms“^ Vs ^ 350 km s“^ (where the 160 km s“^ 
was the maximum shift for a column density of 10^^ cm“^ 
in the frame of the scattering medium). 


® In Figure |Cl| we show that this mechanism produces blue 
bumps for column densities as low as ~ 10^^ cm This 

is because a static shell of gas with 10^^ cm“^ remains 

optically thick to Lya photons emitted close to line centre, and 
can therefore reflect back photons into the expanding shell. 
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Figure 6. Schematic illustration of the origin of photons of the 
blue bump observed in Lya spectra in the shell-model. The out¬ 
flowing shell is surrounded by a static ISM. 


5 DISCUSSION AND CONCLUSIONS 


We have presented a detailed analysis of Fermi acceleration 
of Lya photons across a shock front, a process that was first 
studied analytically by NM88. Our Monte-Carlo radiative 
transfer simulations of the slab model (Figure confirm 
the basic results in NM88. 

In particular, we find that each time a photon crosses 
a shock front that is propagating at Vg, it experiences a 
Doppler boost Vb ~ (0.5 — 0.7)i’s in the gas frame. The 
precise numerical coefficient depends on the opacity of the 
gas to the Lya photons. We also find the blueshift of the 
Lya photons to increase in proportion with the number of 
shock crossings /, i.e. Vb oc (0.5 — L.7)lvs. 

We discussed how our results can help to explain ex¬ 
tended blue wings observed in spectra of radio galaxy 
3C326.1 (which was proposed previously by NM88), but we 
note that our line as a whole appears redshifted by ~ 300 km 
s“^. We reached the same conclusions for the radio galaxy 


B2 0902+34 ( Adams et al.||2009 ). We nevertheless consider 
Fermi acceleration a plausible alternative to the model pro¬ 


posed by Adams et al. (2009), which involves the collapse of 
> 10^^M© of neutral gas. 

Last, we show that Fermi acceleration naturally gives 
rise to blue bumps in Lya spectra, which are difficult to re¬ 
produce with conventional shell-models (at least those with 
reasonable values for the FWHM of the intrinsic Lya line, 
see Hashimoto et al. 2015). We presented a natural extension 


to shell-models in which the shell is expanding into static gas 
that contains a low column density of hQ These models can 
give rise to blue bumps in Lya spectra without abandon¬ 
ing the simplicity of the conventional shell-modej^ In our 


^ If we interpret the shell as dense gas that is swept up by feed¬ 
back processes, then the static gas can be interpreted as gas that 
is being swept up. 

® It is possible that this simply signals the break-down of the 
shell-model, and that radiative transfer through more complex 
gas geometries needs to be studied. There are numerous works 
studying Lya transfer through more complex, simulated gas dis- 


tributions (e.g.|Tasitsiomi, A.|2006||Laursen et al.|2009 

|Barnes et 

|al.|2011||Verhamme et al.|20121|Behrens k Braun|2014| 

Lake et al. 


model, the bump consists of photons that initially streamed 
through the outflowing shell, but which were reflected back 
into this shell by the surrounding static interstellar medium. 
This suggests that blue bumps in Lya spectra are associated 
with outflows that are still confined to the ISM of the galaxy, 
which may represent an earlier stage in the evolution of the 
galaxy. 
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APPENDIX A: ORIGIN OF P(/i) oc ji 

The goal of this appendix is to clarify why the distribution 
of angles at which Lya photons cross the shock front scales 
as P(/i) oc /i (where /i = cos^), rather than isotropically. 

Photons that cross the shock front will on average have 
traversed an optical depth r = 1. Depending on /x, r = 1 
corresponds to a different (frequency dependent) physical 
depth, d{iy), away from the shock front. It is easy to see that 
d oc /i, so that at hxed frequency the volume of gas that 
the photon is likely to have last scattered in is Vcat oc /x. 
If we further assume that the density of scattering events 
is homogeneously distributed in the gas, then this leads di¬ 
rectly to the relation / oc /x, where I is the intensity of 
photons. Therefore P(/x) oc /x, which we showed in Eigurej^ 
and which was found previously by Ahn et al] ( |2001| and 
Garavito-Camargo et al.| ( 2014| ). 

Eigure M shows the average angle at which photons 
cross the shock front, (/x) = dji as a function of 

crossing number 1. As we mentioned in the paper, odd [even] 
number of shock crossings correspond to crossings from left- 
to-right [right-to-left]. This Eigure clearly shows that (/x) ~ 
2/3 for even shock-crossings, while it approaches (/x) ~ 0.5 
for odd-shock crossings when / > 3. 

An explanation for this is offered in Eigure |A2| which 
shows the mean distance from the shock front (scattering 
depth) of the last scattering before crossing the shock as 
a function of the shock crossing number 1. This distance 
follows an inverse exponential distribution. As the photon 
becomes more blueshifted with increasing /, the mean free 
path of the photon also increases, so that for each succes¬ 
sive I the photon is coming from deeper within the slab, 
explaining the increasing trend in Eigure |A2| However, as 
the mean free path increases, the photons which propagate 
the farthest actually propagate all the way through the gas 


1.0 


0.5 


A 

^ 0.0 

V 


-0.5 


I 

Figure Al. Mean projection as a function of shock crossing num¬ 
ber, 1. The solid lines represent the expected values for isotropic 
scattering, while the dashed lines show < /x >= ±2/3. Red trian¬ 
gular points are for odd-l shock crossings, blue square points are 
for even-/ shock crossings. 


^ ^ odd-1 shock crossings 
■ ■ even-1 shock crossings 


slab and exit the simulation. Therefore, the distribution of 
scattering depth is truncated for the next scattering count, 
resulting in a flattening in the growth of scattering depth 
as a function of 1. This is precisely the behaviour we ob¬ 
serve for odd-l in Eigure |A2| The scattering depth distri¬ 
bution is truncated, but the photon blueshift is unaffected. 
Therefore, the gas between the next scattering event and the 
shock front becomes effectively optically thin to the photon, 
and we expect photons to cross the shock front isotropically, 
which is why (/x) approaches 0.5 in Eigure 

Eor even-/ we do not observe this transition to isotropic 
crossing. This is because Ni Nc, so that at a given / and 
blueshift the mean free path of a photon is much lower in 
the right (static) slab and the scattering depth distribution 
for even-/ (i.e. crossings where the previous scattering was in 
the right, static slab) has not yet become truncated. As dis¬ 
cussed above it is the truncation of the scattering depth dis¬ 
tribution that causes the transition to isotropic shock cross¬ 
ing, and so in its absence for even-/ the Lya photons still 
cross the shock front following P(/x) oc /x. 

Einally, we note that the highest shock crossing counts 
that occur in our simulation are rare events, and the re¬ 
sultant poor statistics are responsible for the turnaround 
in Eigure |A2| (recall that we are sampling from an inverse 
exponential distribution). 


Al 


APPENDIX B: INTERSTELLAR DUST 

In dusty media the increase in path length caused by scat¬ 
tering can allow the dust to signihcantly attenuate the Lya 
flux that escapes. Because photons of different frequencies 
have different scattering cross-sections the average number 
of scatterings they undergo will be different. Therefore they 
will have different changes in path-length, and subsequently 
undergo different levels of dust attenuation. 

We take a direct approach to address this issue and 
perform numerical experiments where grey dust is added to 
either the inner, outer, or both shells. The results are shown 
in Figure [BT] It is clear that the differences in the emergent 
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Figure A2. Mean distance between the shock and the location 
of the last scattering before crossing the shock vs shock crossing 
number, 1. d is shown in simulation units, i.e. cells (each slab 
is 128 cells wide). Red triangular points and the red dotted line 
show the odd-1 crossings, while blue square points and the blue 
dot-dash line show the even-1 crossings. 



Figure Cl. Same as Figure 


B1 


but here we vary while 


keeping the other parameters fixed. The prominence of the blue 
peak reduces with HI column density, but a blue bump clearly 
remains even when = 10^^ cm“^. 


APPENDIX C: VARYING THE HI COLUMN 
DENSITY IN THE STATIC SHELL 

Figure [CT] shows that the blue bump remains visible even if 
we further reduce N^i'^ by order(s) of magnitude. 



Figure Bl. Emergent spectrum for a shell-model with 
'i;shell) = (1019 cm-2, IQi^cm-^, 200kms-i) in¬ 
cluding three different dust prescriptions. For each histogram 
'Tdust,i refers to the optical depth of the inner shell, Tdust,o refers 
to the optical depth of the outer shell. The reference dust-free 
spectrum is shown with a heavy black line. 


spectrum caused by the different dust prescriptions are very 
small. 

The addition of dust slightly enhances the blue bump. 
This occurs because once a photon is Fermi accelerated the 
optical depth of the HI shells to the photon is vastly re¬ 
duced and the photon easily escapes, avoiding or reducing 
the number of further scatterings. Thus, a Fermi acclerated 
photon has a shorter total path-length so that the presence 
of dust affects it less than an unaccelerated photon. This 
results in a slightly enhanced blue peak in the normalised 
spectra shown in Figure [BT] 
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